Optical selection rules for electron-hole pair excitation in 3D topological insulators 
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Experiments using ARPES, which is based on the photoelectric effect, have shown that the surface 
states in 3D topological insulators (TI) are helical. Here we consider Weyl interface fermions due to 
band inversion in narrow-bandgap semiconductors, such as Pbi-ccSn^Te and Bii-^Sbx. We deter- 
mine the optical selection rules of electron-hole pair excitation by means of the solutions of the 3D 
Dirac equation. While electron- hole pairs in graphene are generated through intraband transitions, 
we show that in 3D TI they are generated through both intraband and interband transitions. For 
their analysis, we calculate explicitly the electric dipole matrix elements by means of bandstruc- 
ture calculations for Pbi-ccSn^Te. While in graphene the pseudospin helicity operator is given by 
h g — — (1/ cr • p, we define the spin helicity operator in 3D TI as /iti = — (1/ (cr x p) ■ z, 
where z points perpendicular to the interface. Our results are crucial for future opto-electronic 
devices based on 3D TI. 
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The 3D topological insulator (TI) is a new state of mat- 
ter with unusual properties that manifest themselves on 
the surface or at the interface of narrow band gap ma- 
terials where topologically protected gapless (metallic) 
surface/interface states appear within the bulk insulat- 
ing gap [T]-[5]. These surface/interface states are charac- 
terized by the linear excitation energy of massless Weyl 
fermions. Protected by time-reversal symmetry, the spins 
of the Kramers partners are locked at a right angle to 
their momenta due to spin-orbit coupling, making them 
robust against perturbation and scattering [IJ |2j |6j |7]- 
This means the larger the spin-orbit coupling, the better 
is the protection. 

TIs can be found in a subclass of narrow-gap semicon- 
ductors where the band gaps are smaller than the typical 
energy scale of relativistic corrections to the bandstruc- 
tures. The 2D TI exhibiting the quantum spin Hall effect 
was experimentally observed in HgTe quantum wells [8], 
where the interface of the well lead to a sign change of 
the mass term in the 2D Dirac equation, thereby creat- 
ing a pair of gapless helical edge states inside the bulk 
energy gap. The origin of this hallmark for both 2D and 
3D TIs is the Rashba spin-orbit coupling [9 J . Helical edge 
modes have also been reported in a heterostructures of 
InAs/GaSb [TO] . Many materials have been proposed as 
3D TIs [llj. Because of the presence of a single Dirac 
cone with fixed spin direction at the surface, the main 
feature of strong TIs [12, 13J, the materials Bi2Se3 and 
Bi 2 Te3 are currently being widely studied [5, 14J. 

The heterostructures of compound semiconductors 
such as Bii-^Sb^ and Pbi-^Sn^Te exhibit a strong topo- 
logical phase |7J. In Bii-^Sb^, with increasing x the gap 
at the L point decreases and the L + and L~ bands cross 
at x = 0.04. After further increase in x, the gap reopens 
with inverted band ordering. In recent Angle Resolved 
Photoemission Spectroscopy (ARPES) experiments, sur- 
face states with an odd number of crossings (five surface 
states crossings) between Y and M at the Fermi energy 



have been observed in this material with x = 0.1 |10| . 
The pure PbTe has inverted bands at the band gap ex- 
trema with respect to SnTe. In Pbi-^Sn^Te, initially 
increasing the concentration of Sn leads to a decreasing 
band gap. At around x = 0.35 (at Helium tempera- 
ture), the bands cross and the gap reopens for x > 0.35 
with even parity L + band and odd parity L~ band being 
inverted with respect to each other [15J. The band inver- 
sion between PbTe and SnTe results in the existence of 
helical interface states fT6HT8] . which according to Volkov 
and Pankratov can be described by a 2D mass-less Dirac 
equation, i.e. the Weyl equation |19| . 

In this Letter we investigate the optical selection rules 
governing the low energy excitation of electron-hole pairs 
around a Dirac point in a 3D TI. Note that the selection 
rules obtained here are different from the selection rules 
in ARPES experiments, which record the number of pho- 
toelectrons as a function of kinetic energy and emission 
angle with respect to the sample surface. A number of ex- 
periments have shown the existence of the helical surface 
states in 3D TI However, a clear understanding of 

helical surface/interface states in terms of solutions of the 
Dirac equation is still lacking. Here we show that the heli- 
cal states in 3D TI naturally emerge from the solutions of 
the 3D Dirac equation in the presence of a band inversion. 
As an example, we consider the alloy Pbi-^Sn^Te, which 
has topologically nontrivial interface states under appro- 
priate doping level. Our results are valid for all strong 3D 
TIs. Pbi-^Sn^Te has a rocksalt type crystal lattice with 
four non-equivalent L points located in the center of the 
hexagonal facets on [111] axis. The valence and conduc- 
tion band edges are derived from the hybridized p-type 
orbitals at the L point [20J. Its end species have inverted 
band character, L + character of PbTe band switches to 
L~ character of SnTe band and vice versa as shown in 
Fig. [T] We choose the z-axis to point in direction of the 
gradient of the concentration Vx. At the two band ex- 
trema, the low energy Hamiltonian is described by a 3D 
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Figure 1: a, Band inversion in between two end members 
in Pbi-^Sn^Te. b, Energy spectrum of the inverted contact. 
The solid lines are Weyl states and dashed lines are additional 
branches appear for contact thickness I > l Q . At the band 
edges, dispersion is highly linear. 

relativistic Dirac equation whose solutions are localized 
near the z = plane where the band crossing occurs, 
which defines the interface. Dispersion is nearly linear 
owing to the large band velocities of v± = 8 x 10 5 m/s 
and vn = 2.24 x 10 5 m/s with a small gap [19J. Such 
peculiar properties at the interface result in a small lo- 
calization length l Q of the wave function of interface states 
along the z-axis. Due to the absence of a center of inver- 
sion, a Rashba-type spin-orbit coupling is present, which 
is automatically taken into account in the framework of 
the Dirac equation. We also present the details of our ab- 
initio calculation of the bandstructures of these materials 
in the supercell Brillouin zone obtained by doubling the 
lattice parameters in each direction. Analysis of the al- 
loy band structures is usually complicated due to folding 
of the bands from neighboring Brillouin zones, making it 
difficult to map the calculated bandstructures onto the 
bandstructures obtained from momentum-resolving ex- 
periments. The analysis is further complicated by the 
presence of impurity bands inside the normal bulk en- 
ergy gap. The interface states sometimes overlap with 
bulk energy states. Therefore, we unfold the band struc- 
tures along the [111] direction in order to shift the band 
crossing from the T point, as seen in the supercell Bril- 
louin zone, to the L point in the primitive cell Brillouin 
zone pTll22]. 

The energy spectrum of Pbi-^ Sn^Te near the 
band crossing is described by a Dirac Hamiltonian. 
Within the envelope function approximation, neglecting 
the far band terms, the interface is described by the two- 
band Hamiltonian 

H -{a.p -A(z) ) W 

where a are the Pauli matrices, p is the momentum oper- 
ator with projections — zfti^V^and —ihvnVn on the xy- 
plane and along the z axis, respectively, A(z) = e g (z) /2 
is the gap energy parameter that varies along the z- 
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Figure 2: Brillouin zone for rocksalt type crystal with space 
group Fm3m. There are four inequivalent L points at the 
center of the faces on the surface of the Brillouin zone. The 
growth direction is along [111] and is chosen to point along 
the z-axis . With the appropriate level of doping by Sn atoms 
in PbTe, band gap goes to zero at L point with a linear exci- 
tation energy that traces out a cone in the 2D Fermi surface 
parallel to the face of Brillouin zone that is perpendicular to 
the growth direction. 

axis. Note that A (z) breaks the translational invari- 
ance of the crystal in the z direction and has a symmetry 
A (z) = — A (— z). The inhomogeneous structure is syn- 
thesized by changing the composition along one of the 
[111] axes, which can be chosen along the z direction. 
After a unitary transformation of the Hamiltonian, the 
time-independent Dirac equation can be written as 

/ iA (z) + a • p \ / x- \ = 

V ~iA (z) + a • p J \ X + J 

(-*■<•»(£) P) 

where Xt are ^ ne ^ wo component spinors. The potential 
cp (z) (work function) describes the variation of the gap 
center. For simplicity we consider the case cp(z) = 0. 
From equation the two component spinor x± satisfies 

(p 2 + u±(z,a z )-s 2 )x±=0 (3) 

where u±(z,a z ) = A 2 ± hv\\cr z ^. In its origin, the 
linear Weyl spectrum (k±) = ±hv_\_k_\_&t k± = is 
approximately equal to the soliton spectrum in the ID 
Peierl's insulator. This implies that A (z) can be cho- 
sen to be A(z) = A (oo) tanh (z/l). Interface states are 
localized along the z-axis with the localization length 
l Q = hv_\_/A(oo). For l Q < I, additional branches 
with finite mass appear. There are several solutions at 
e 2 > A 2 (oo) which are localized at the contact. For 
l > Z, only Weyl solutions exist. We focus on the case 
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when l Q > I. Then we have only zero-energy solutions, 
which correspond to the Weyl states and are given by 



* ± = c 
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where C is a normalization constant, = (k x , k y , 0) 
and e^ l ° = kx J^ ky . These solutions have eigenenergies 
(^-l) = For A (z) to vanish at the inverted 

contact, it can be seen from Equation Q that x+ an d 
X- must have only non-zero spin down and spin up com- 
ponents, respectively. Each spinor is a mix of spin states 
from the L~ and L + bands. Each spinor at L T band 
can be represented with the spin up states from the L~ 
band and spin down states from the L + band for both 
the positive and the negative energies. The motion of 
the particle at the inverted contact is separated into free 
motion in (p xi p y ) and confinement in p z . A remarkable 
property of equation ([2| is the presence of the zero mode 
(Weyl mode) localized around z = 0. It is this mode 
that has a locked spin structure. In order to understand 
the direction in which the 4-spinors point, we have to 
transform the solutions back to the original basis of the 
Hamiltonian in Equation ([!]). With the transformation 
ty± = U^± where U is the unitary matrix used for the 
basis transformation, the Weyl solutions are 
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where the global phase factors of 



(5) 

are ignored. As 



a matter of fact, the solutions are helical. In the rep- 
( \ 

resentation \£± = ( f f+ ] , the two spinors and 




ip± refer to the L~ band and L + band respectively. It 
is important to note that the solutions for each sign of 
the energy are mixed helical states from both bands. 

In this representation, the spin directions reveal 
themselves clearly: Defining F(r) as the position- 
dependent part, the spins of the 

• (0±7T/2) \ 

F(r) and = 



2-spinors ip± 



i (g£[/2) 



(e±7T/2) 



F(r) 



point perpendicular to k^ owing to the shifts. This 

important property will be emphasized with the defini- 
tion of a new helicity operator for 3D TI later in this 
Letter. For the states corresponding to the positive eigen- 
value at both bands L T , the directions of the spin vectors 
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Figure 3: Effect of Rashba spin-orbit coupling. Spin vector 
in the L~ band and spin vector in the L~ band 

are perpendicular to both the z-axis and p ± . Orientations 
of the up spin component and down spin component are in 
opposite directions to each other, a and b are for the electron 
with energy, e + = +hv±k± and c and d are for the hole with 
negative energy, e~ — —hv±k±. 



where the signs =p refer to the L T bands, associ- 
ated with the two component spinors point oppositely to 
each other. For an asymmetric scalar potential V ap- 
plied to a semiconductor heterostructure, the inversion 
symmetry is broken, which leads to the Rashba spin- 
orbit coupling [9j [23]. Here in the case of the interface 
of a 3D TI we have an asymmtric potential V T = ±A, 
whose sign depends on the band L T . This results in a 
band-dependent Rashba spin-orbit coupling. For the pos- 
itive solution (electrons) the Rashba spin-orbit coupling 
has the form H R = -X R cr • (p x VV T ) = -A^VV^ • 
(cr x p) = TA^VA • (cr x p) = T^R^t ( a *Py ~ a yP*)i 
where Xr > is the Rashba spin-orbit coupling constant. 
Similarly, for the negative solution (holes), the Rashba 
spin-orbit coupling is given by Hr = \r<T'(p x V7 t ) = 
±Xr^ (cr x p y — CFyPx). In both cases each spin is 
perpendicular both to the momentum and to the poten- 
tial gradient direction, i.e. the z-axis (see Fig. |3|. Our 
analytical findings are consistent with the spin density 
functional calculation (DFT) done by Zhao et al. |24| . 

In a strong TI the surface Fermi circle encloses an odd 
number of Kramers degenerate Dirac points. The sim- 
plest case is a single Dirac point. The surface electronic 
structure of a 3D TI is similar to a graphene but here we 
have a single Dirac point rather than having four Dirac 
points (2 valley x 2 spin) |7J. In order to determine the 
Kramers partners explicitly, we rotate the phase of each 
of the two component spinors by an angle 7r in the 2D 
interface plane (see Fig. [4|. This rotation yields the 
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Figure 4: Kramers partners obtained by a n rotation of each 
two component spinor. Spin direction and momentum di- 
rection of each partner are nipped with respect to the spin 
direction and momentum direction of their respective part- 
ners in original basis, a and b are for positive eigen value and 
c and d are for negative eigen value. 



Kramers partners ip± (0 - 



and (0 + tt) 
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.(e^Tr/2) 



F(r) 



F(r), where global 



phase factors of e ±2 ^ have been ignored. Both the spin 
and momentum direction of each partner component is 
flipped by an angle tt (Fig. |4|. This provides a theoreti- 
cal hallmark of the existence of Kramers partners in 3D 
TL 

The relevant quantity to characterize the eigenfunc- 
tions is their helicity defined as the projection of the mo- 
mentum operator defined along the spin direction. Heli- 
cal properties of solutions given by equation ([5| are com- 
pletely general to the 3D TL Our analytical results show 
that the spins point in perpendicular direction to the 
momentum. In the case of free neutrinos in 3D space, 
the standard helicity operator h n = — (1/ a • p for 
the spin S = ha/2 can be used. Similarly, in the case 
of graphene the helicity for the pseudospin is given by 
h g = — (1/ |rl|) <t • p. However, in the case of 3D TI 
this definition is not useful, because the spin points per- 
pendicular to the momentum. Therefore, since we know 
that the Rashba spin-orbit coupling is responsible for the 
helicity in 3D TIs, we define the 3D TI helicity operator 
as hn = - (1/ \p ± \) (cr ± x p ± ) ■ z, where a± = {a x ,a y } 
is the two dimensional vector of Pauli matrices in the 
xy-plane. h^i commutes with the effective single-band 
Hamiltonian H for the L~ or L + band if the expecta- 
tion value of the spin along the z-axis is zero, which is 




Figure 5: Dirac cones at the crossing between the L and 
L + bands. The orange (blue) arrows denote interband (in- 
traband) transitions. In both types of transitions the spin is 
conserved. Spins are always tangential to the circle of radius 
p±on the Fermi surface. The color of the cones represents the 
helicity of the Weyl fermions in the band. In the L~ band 
the Weyl fermion with positive (negative) energy has positive 
(negative) helicity, which is denoted by a green (blue) color 
cone. The helicities are inverted in L + band. Note that both 
cones belong to the same time reversal invariant momentum, 
in contrast to the K and K' points in graphene. The thresh- 
old energy above which an electron-hole pair can be excited is 
given by the difference in energy between the Dirac point and 
the Fermi level. The reason for this threshold is the Pauli ex- 
clusion principle, which does not allow for double occupation. 



valid in our case because the spins lie in the x?/-plane. 
The eigenvalue of the operator h^i for the positive en- 
ergy is (+1) for the L~ band and (—1) for the L + band, 
i.e. /iti'0+ =F = (±1)^+^ The eigenvalues of h^\ are re- 
versed in the case of negative energy, i.e. (—1) for the 
L~ band and (+1) for the L+ band, h Tl ^ = (=pl) ^ 
(see Fig. |5|. 

Now we proceed to calculate the optical selection rules 
for the excitation of electron-hole pairs, keeping in mind 
that the Dirac equation provides an effective description 
of the two-band system consisting of the L T bands. It 
is interesting to compare our results obtained for 3D TI 
with graphene. Let us represent the four spinor s \£± as a 
linear combination of the two component spin states from 
L T bands including the periodic part of the Bloch func- 
tion from each band, \£± = a±u% F (r) + P±u^ F (r), 

/ ±e ~i^ 1 \ ( c -i^l \ 

where a± = I i (e± 2L /2 1 I and (3± I , (9Tff/2 ) I . 

The Bloch functions in the two materials on either side 
of the interface must be similar in order for the effective 
mass approximation to be valid. An obvious condition 
is that they must belong to the same point in e n (k). 
We calculate the low energy transitions around an L 
valley that is lifted up along the z-direction from the 
other three L valleys. It is important to note that in the 
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material we are investigating, with the proper choice of 
uniform strain, composition and layer width, there exist 
practically gapless helical states for the [111] valley. In 
unstrained Pbi-^Sn^Te, band inversion occurs simulta- 
neously at four L points and the phase is topologically 
trivial. For most experiments, in a structure with inter- 
face thickness d « 10 nm, dispersion of the [111] valley 
states can be assumed to be gapless while the states in 
the oblique valleys are gaped [25J. This realization is 
analogous to the situation in 3D topologically nontrivial 
materials such as Be2Tes/Bi2Se3 where the difference in 
gaps is expected. The helical properties of the localized 
spin states of 2D Weyl fermions in 3D TI are different 
from the ones in graphene where the solutions are eigen- 
functions of the tight-binding Hamiltonian. Dirac cones 
at the L~ band and L + here resemble the Dirac cones at 
the K and K points in graphene; however, both L points 
in a 3D TI belong to the same k-point in the Brillouin 
zone. Under the dipole approximation, the transition 
matrix elements in a 3D TI are given by 



\e -p\ *_) = -4ft sin 



4zftcos 



e • z 



ie k\ 



)^e-Glm(al + a L -) (6) 

G 



where e is the direction of the polarization, G is the recip- 
rocal lattice vector and a L ± are the expansion coefficients 
of the periodic part of Bloch's functions. The envelope 
functions are properly normalized. The transition prob- 
ability can be easily calculated using Fermi's golden rule. 
Equation (|6| has a number of interesting features. When 
we carefully analyze right hand side of the equation (J6|, 
it is evident that the first term corresponds to transi- 
tions within the L~ band and within the L + band (i.e. 
intraband) and the second term corresponds to transi- 
tions between L~ band and L + band (i.e. interband), 
which are shown in Fig. [5] For the parallel polarizations 
of light, intraband as well as interband transitions occur 
simultaneously. The same is true for the perpendicular 
polarization of light with different strengths of the transi- 
tion. It is also interesting to note the features associated 
with equation ^ as a function of polarization angles. On 
the one hand, for — = 0, transitions are purely inter- 
band determined by the crystal symmetry. On the other 
hand, for — 6 = n, transitions are purely intraband 
with a contribution coming from the envelope function. 

The relative weight of the transitions in equation ([6| 
requires a complete calculation of the bandstructures. 
Therefore, let us calculate the bandstructures and in- 
volved dipole matrix elements. Fig. [6] shows our calcula- 
tions of the complete bulk bandstructures of Pbi-^Sn^Te 
at 37.5% doping by Sn impurities in a supercell Brillouin 
zone using density functional theory within PAW approx- 
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Figure 6: Bulk bandstructure of Pbi_ x Sn x Te at x = 0.375 
doping level including spin-orbit coupling. In the experiment, 
the crossing has been reported nearly around x = 0.35 at He- 
lium temperature [15] . In the supercell Brillouin zone bands 
are folded from the neighboring Brillouin zones into the first 
Brillouin zone (a and b). A minimum band gap of 1.8 meV 
appears at the F point of the supercell Brillouin zone. In fact, 
the minimum band gap occurs at the L point as in the band- 
structures of its parent compounds PbTe and SnTe. This is 
seen as the bands are unfolded in the primitive cell Brillouin 
zone, as shown in Figure c, where the bands are unfolded 
from the F to the L point in A direction. Lsc and Lpc are 
L points in the supercell Brillouin zone and the primitive cell 
Brillouin zone, respectively. The solid green color in Figure 
c denote the spectral functions calculated using the recipe 
described in reference [21] . Bands of opposite parity nearly 
cross at around 67 meV below the Fermi level at the L point 
where we observe Dirac cones, Figure c. The unfolded band- 
structure facilitates the accurate measurement of the band 
crossing. A threshold energy of around 67 meV or above is 
required to excite an electron-hole pair at the L point (d) at 
37.5% doping level. 



imation as implemented in VASP (Vienna ab initio Sim- 
ulation Package) [26-28J. We unfold the bandstructures 
along the F point to the L point using the unfolding 
recipes developed by Popescu and Zunger [2T] , 

There is a small band gap of 1.8 meV at the F point. 
In the unfolded structure the same amount of band gap 
appears at the L point as in folded structure. This is 
required by the conservation of energy. In the unfolded 
structure the bands around the L point are almost linear, 
which is best described by Weyl fermions. Two bands 
with opposite parity nearly cross each other at around 
67 meV below the Fermi level at the L point. In a recent 
experiment on Bi2Se3 film grown on sapphire, the Dirac 
point was observed at 135 meV below the Fermi level 
[29]. In the ARPES measurement for the verification of 
PbBi2Te4 as 3D TI, the Dirac point was observed at 470 
meV below the Fermi level [30J. In the unfolded structure, 
the valence band maximum is derived from the p orbitals 
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of Pb and Sn hybridized with the s orbital of Te and 
conduction band minimum is derived from the s orbitals 
of Pb and Sn hybridized with the p orbital of Te. The 
valence band maximum and conduction band minimum 
have opposite parity, thus making interband transitions 
allowed. The interband dipole matrix elements at 37.5% 
doping are calculated to be 3.82 x 10 -26 Js/m for the 
perpendicular polarization and 1.12 x 10 -25 Js/m for the 
parallel polarization of the light. The localization length 
l Q for the Dirac Fermion along z-axis can be obtained 
using the calculated band gap of 350 meV including spin- 
orbit coupling for PbTe (Experimental band gap is 190 
meV). Using the band velocity, v\\ = 2.24 x 10 5 m/s, 
l Q is found to be 0.84 nm. This length measures the 
characteristic scale of confinement for the Weyl states 
that exists at the interface. 

The dipole matrix elements given by equation ([6| are 
used to calculate the transition probability given by the 
Fermi's golden rule, shown in Fig. [7] For the parallel 
polarization of the light, the intraband dipole matrix ele- 
ments are a little larger than the interband dipole matrix 
elements, as can be seen in Fig. [7^i. For the perpendicular 
polarization of the light, intraband transitions are dom- 
inant, as can be seen in Fig. [TJd. For intraband transi- 
tions, the transition probability is maximum when the 
momentum in L~ and in L + bands are opposite to each 
other. 3D polar plots are shown in Figs. ^ and d. They 
visualize the transitions as a function of polarization an- 
gles in k± = k Q (x cos 6k + y sin 0k). The inital and final 
state have momentum angles 0^ = and k = , re- 
spectively. For = and 6 =0, transitions are purely 
interband occurring from L~ to L + and vice versa with 
slightly higher probability in xz-plane. This type of tran- 
sition is absent in graphene. For 6 = and 6 = 7r, tran- 
sitions are purely intraband occurring only in xz-plane, 
transitions along 7/-axis are completely suppressed. This 
is because the momentum of the initial state is aligned 
along the x-axis. The transitions occurring along z- axis 
are governed by the envelope function. 

The topological interface states occurring in 3D topo- 
logically nontrivial phase are probed using experimental 
techniques such as ARPES and scanning tunneling mi- 
croscope [3Tj . The states are shown to be helical. In 
this Letter, we show that these states can be derived di- 
rectly by solving the Dirac equation in 3D. The nature of 
the helicity of the states is revealed in the solutions. We 
propose a new helicity operator that identifies the helicity 
of Weyl fermions in general in all 3D TI. Our theoretical 
proof of the existence of Kramers partners in 3D TI and 
helical interface states is very important for the further 
study of the non-trivial topology of the electronic prop- 
erties. From our derivations, it has become clear that 2D 
helical interface states occurring at the interface of 3D TI 
are different from the helical states occurring in graphene 
in a number of ways. Most importantly, the optical selec- 
tion rules proposed for the Weyl fermions in this Letter 
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Figure 7: Transition probability calculated using the Fermi's 
golden rule using the dipole matrix elements given by equa- 
tion |6|. Units are arbitrary. For the incident light polarized 
parallel to z-axis, both the intraband and interband transi- 
tions are significant (a) but the strength of the intraband 
transitions is higher. For the incident light polarized perpen- 
dicular to the z-axis, the transitions are mostly intraband (b). 
3D polar plots are shown (c and d) to visualize the transition 
probability occurring on a sphere as a function of polarization 
angles. For = and =0, transitions are purely interband 
and for = and = tt transitions are purely intraband. 



show that the transitions are both intraband and inter- 
band at the L point, in contrast to graphene where the 
transitions are purely intraband (Fig. |5j. These selection 
rules (equation (J6|) we found are of great importance for 
the excitation of electron-hole pairs in low power opto- 
electronic applications. 
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